EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



We search for lepton flavour violating events (e/i , er and fir) that could be directly 
produced in e + e~ annihilations, using the full available data sample collected with 
the OPAL detector at centre-of-mass energies between 189 GeV and 209 GeV. In 
general, the Standard Model expectations describe the data well for all the channels 
and at each y/s. A single e/i event is observed where according to our Monte Carlo 
simulations only 0.019 events are expected from Standard Model processes. We obtain 
the first limits on the cross-sections cr(e + e~ — > e/i, er and /rr) as a function of y/s at 
LEP2 energies. 
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1 Introduction 

Within the minimal Standard Model (SM), the fermion mass matrices and the mech- 
anism of electroweak symmetry breaking remain unexplained. The conservation of 
lepton number separately for each generation has no strong theoretical basis. In 
addition, recent data JT| present evidence for neutrino oscillations which necessar- 
ily violate lepton-flavour symmetry. Beyond the SM, lepton flavour violation (LFV) 
can occur in many supersymmetric (SUSY) extensions. An example is the SO(10) 
SUSY GUT model where both the left- and right-handed supersymmetric lep- 
ton partners induce LFV, but none of the existing models predicts a measurable 
effect at LEP2 energies. Experimentally, no evidence for direct LFV has been re- 
ported so far. Upper bounds for muon decays are BR(jU~ — > e~7) < 1.2 x 1CT 11 
and BR(jU + — > e + e + e~) < 1CT 12 ||. Searches for neutrinoless r decays such as 
r + — > e + e + e~ and r + — > yield upper limits BR(Z — > er) < 5.4 x 10 -5 and 

BR(Z — > fir) < 7.1 x 1CU 5 at 90% confidence level (CL). Direct searches in e + e~ 
annihilations at the Z peak performed by the LEP experiments || yielded 95% CL 
limits BR(Z — > e/i, er, /it) < O (1) x 1CT 5 . In this paper we search for the e/x , er 
and \it final states in e + e~ collisions at centre-of-mass energies between 189 GeV and 
209 GeV. 



2 Data sample and event simulation 

The OPAL detector is described in detail in ||. In the present analysis, the silicon 
micro vertex detector, the central tracking chambers, the electromagnetic calorimeter, 
the hadron calorimeter and the muon chambers were required to be fully operational. 
The full data sample collected since 1998 at = 189 GeV and above is analysed. 
The corresponding collected integrated luminosities are shown in Table 1. 

Track reconstruction is performed by combining the information from the silicon 
micro vertex detector, the vertex drift chamber, the large volume jet drift chamber 



3 



and an outer layer of drift chambers for the measurement of the z coordinate^ The 
OPAL electromagnetic calorimeter (ECAL) consists of a barrel part covering the 
region | cos^l < 0.82 and two endcaps covering the region 0.82 < | cos#| < 0.98. A 
set of forward detectors provide complementary coverage for 9 > 25 mrad. 

Detector acceptance and reconstruction efficiencies for processes under study are 
evaluated with two different methods based on Monte Carlo simulations. They are 
cross-checked with a third method that uses data only, where possible. The first 
method consists of generating e/i , er and \xr signal events with the EXOTIC [[7] 
generator using isotropic angular distributions. The second method uses Standard 
Model lepton pair final state events. The lepton pairs were simulated using kk2f [El 



and KORALZ |K| for rrfr) and (ifi{j) and using BHWIDE |I| and TEEGG for 
ee(7). Events passing the preselection are mixed to obtain e/i , er and fir topologies. 
The mixing assumes a uniform detector response in phi, and events are mixed only if 
the polar angles of their thrust axes are within 2 degrees of each other. The mixing 
is performed by rotating the momentum components in the plane transverse to the 
beam axis of a lepton from one event to match the lepton replaced in the mixed 
event. For each channel, the efficiency is estimated in bins of cos# and then averaged 
giving the same weight to each bin, in order to obtain a value that corresponds to a 
uniform angular distribution of the particles in the final state. The third method is 
the same as the second method but it uses a high purity sample of lepton pair events 
from the real data. The efficiencies after the final event selection (described in the 
next section) are almost independent of the centre-of-mass energy. They are shown 
in Table 1, after averaging for different centre-of-mass energies, together with their 
statistical and systematic uncertainties added in quadrature. 

The SM background contributions are evaluated with large samples of events 
processed through a full simulation of the OPAL detector || and analysed using 
the same reconstruction and selection programs as applied to the data. The SM 
sample comprises lepton pair final states with initial and final state radiation, qq(7) 
events generated with PYTHIA [Tj|, a full set of four-fermion final states generated 



with grc4f and KORALW |15| and gamma-gamma scattering events generated 



using HERWIG |16| and PHOJET [17]]. SM Monte Carlo processes were generated at 
y/s = 189 GeV, 192 GeV, 196 GeV and in steps of 2 GeV from 200 GeV to 208 GeV. 
The total generated SM sample corresponds to more than 500 times the integrated 
luminosity of the recorded data for r + r~ , /i + /i~ , qq and 4-fermion final states and 
to about 50 times the integrated luminosity of the recorded data for e + e~ final states. 



3 Event selection 



Events with at least two and at most 8 measured tracks and no isolated photons in 
the electromagnetic calorimeters are considered for the analysis. Isolated photons 

1 The OPAL coordinate system is denned so that the z axis is in the direction of the electron 
beam, the x axis points towards the centre of the LEP ring, and 9 and <p are the polar and azimuthal 
angles, defined relative to the +z- and + re-axes, respectively. In cylindrical polar coordinates, the 
radial coordinate is denoted r. 
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are defined as in [18]. Each track should be consistent with originating from the 
interaction point with a measured momentum, P tr k, exceeding 250 MeV. An energy 



flow algorithm |19], E0[ is used to measure each track and cluster energy, and to correct 



for possible double counting. The cone jet finder ||21|| , with a cone half angle of 15° 
and a minimum cone energy of 10% of the beam energy, Eb ea m, is used and 2 jets are 
required in each event. Here, a single isolated track can form a jet if it has more energy 
than the minimum required energy in the cone. To gain in the energy resolution and 
remove potential background, namely Bhabha and gamma-gamma scattering, the 
momentum vector of each jet must satisfy | cos#| < 0.82. 

All three search channels have in common the feature that one of the event hemi- 
spheres should consist of a single electron or a single muon with a measured momen- 
tum close to Ebeam- Events are selected for further analysis if: 



• The total measured energy outside the two cones defining the jets is smaller 
than 0.10 x Eb ea m, the sum of the two jet energies is larger than Eb ea m and the 
event thrust is greater than 0.95. To have sensitivity to events produced at 
effective centre-of-mass energies lower than the actual a/s, we accept events for 
which the total measured energy in both ends of the forward detectors is up to 
10 GeV. 

• One of the two jets has its energy greater than 0.8 x Eb eam - That jet should 
contain a single isolated track, defined as a track with a momentum greater 
than 0.5 x Eb eam located in a cone with a 10° half opening angle. The 15° 
cone should not contain any other track that has a momentum in excess of 
0.03 x E beam . 



After these cuts most of the hadronic final states and the 77 scattering events are 
rejected and the remaining event sample consists of 92% wide-angle Bhabha scattering 
events and ~ 8% and r + r~ events. The distributions of the visible energy 

(defined as the sum of all measured track energies), the total measured energy in 
the electromagnetic calorimeter and the cos# of the total momentum are shown in 
Figure |I[ Here, the total momentum is defined as the vector sum of all measured 
momenta. The SM expectation describes the measured data well. 

A track is considered to be the electron candidate if it has an associated en- 
ergy in the electromagnetic calorimeter, E EC al, within 20% of Eb eam and if the ratio 
EECAL/Ptrk is greater than 0.7. The track should also have a characteristic ionisation 
in the tracking chambers consistent with an electron hypothesis. A track is considered 
to be a muon candidate if it has matching hadron calorimeter and muon chamber hits 
and if the ratio E EC AL/Ptrk is less than 0.1. The following cuts, common to all three 
search channels, are applied: 

1. To reject the Bhabha scattering, events are selected if the total measured energy 
in the electromagnetic calorimeter is less than 1.6 x Eb eam - 
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Figure 1: The measured visible energy, the total measured energy in the electro- 
magnetic calorimeter and the cos# of the total momentum for all data (points) as 
compared to the Standard Model Monte Carlo expectation (solid lines). The small 
asymmetry in the total momentum distribution is caused by the fact that the inter- 
action vertex is not the geometrical center of the OPAL detector. They are about 
~ 1 cm away from each other. But this effect is well modelled in the Monte Carlo. 
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2. Events that have a reconstructed electron candidate in the transition region 
between the barrel and endcaps (0.75 < |cos#| < 0.85) are rejected, ensuring 
that the remaining electron candidates have a good energy resolution. 

3. To reduce the e + e~ and /x + /x~ background further, events that have two muon 
or two electron candidates are rejected. This reduces the er and /xr efficiencies 
by about f 5% each, rejecting events with r — > uue or uuji where the final state 
electron or muon satisfies the requirement of an isolated track associated with 
an electron or a muon. 

4. To suppress the r + r~ final states, events with two isolated tracks are rejected 
if each track has a measured momentum less than 0.30 x E bcam . They are also 
rejected if the opening angle between the two isolated tracks is less than 160 
degrees. 

The distribution of EECAL/Ebeam for electron candidates and the distribution of 
Ptrk/Ebcam for muon candidates are shown in Figure 2. These distributions are ob- 
tained for the selected events after cut 1 above and show good agreement between 
data and SM expectations. Using Monte Carlo lepton pair final states, we find that 
the resolution on Eecal is 3.5% for electrons and the resolution on P tr k is 11% for 
muons. 

Events that survive the above cuts are subject to specific e/x , er and /xr selection 
cuts. These cuts are optimised to minimise the dependence on y/s and to reject SM 
lepton pair final states while keeping reasonable efficiency for e/x , er and /xr . 

• Selection for the e/i. channel: The event is required to have only two isolated 
tracks with opposite charge where each of the tracks belongs to a different jet. 
One of the tracks should be identified as an electron and its measured energy 
should be within 15% of the beam energy. The second track should be identified 
as a muon candidate with a measured momentum within 20% of the beam 
energy. The total event momentum is required to be less than 0.25 x Eb ea m- 
In the case where the opening angle between the two tracks is less than 170°, 
we require in addition that the missing momentum should be pointing to the 
forward part of the detector (|cos#| > 0.9). Missing momentum caused by 
undetected neutrinos from a r decay would, for the selected events, point to the 
barrel part of the detector. 

• Selection for the er channel: Here we apply tighter cuts on the electron 
candidate. The electron should be identified as an isolated track with |E EC al — 
Ebeam | /Ebeam < 0.10. The energy of the recoiling jet against the electron is 
required to be less than 0.75 x Eb ea m, and the total event momentum should be 
larger than 5 GeV and point to the barrel region of the detector. 

• Selection for the \it channel: Here we require an identified muon as an 
isolated track with |P tr k — Ebeam | /Ebeam < 0.15. The energy of the jet recoiling 
against the muon is required to be less than 0.75 x Ebcam, and the total event 
momentum should be larger than 5 GeV and point to the barrel region of the 
detector. 
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Figure 2: The E EC AL/Ebeam distribution for electron candidates and the Ptrk/Ebeam 
for muon candidates for all data (points) as compared to the Standard Model Monte 
Carlo expectation (solid lines). 
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Table 1: Integrated luminosity and efficiency as a function of y/s. 



v/i(GeV) 


Lumi. [pb 1 




e/i [%} 


er 


%] 




fir [%} 


189 


174.6 


56.3 ±0.7 


24.6 ±0.5 


22.2 ±0.6 


192 < v/i < 200 


103.6 


56.0 ±0.7 


24.2 ±0.5 


21.8 ±0.6 


200 < y/s < 209 


322.2 


55.8 ±0.7 


23.9 ±0.5 


21.8 ±0.6 



Table 2: Selected data events versus SM expectations. "No— £ + £~" stands for the cut 
against events having a pair of electron (e + e~) or a pair of muon (n + n~) candidates 
in the final state. 





e + e 






Other 


Total Background 


Data 


Selected 


20745 


1359 


564 


24 


22683 


23164 


I]E E cal < 1-6 E beam 


2275 


1359 


520 


23 


4185 


4201 


No - t+fr 


559 


57 


67 


21 


704 


713 


efjL Candidates 








0.015 


0.004 


0.019 


1 


er Candidates 


4.010 


0.017 


0.520 


0.004 


5.01 


5 


/j,t Candidates 


0.017 


5.901 


8.400 





14.3 
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With our lepton identification criteria, a hadron from a single-prong tau decay has 
a 3% probability to be misidentified as an electron and a 1.6% probability to be 
misidentified as muon. 



4 Results 

After the normalisation of the Monte Carlo backgrounds to the integrated luminosity 
of the data, the effects of the selection criteria are summarised in Table 2. The 
contributions to the SM background from processes other than lepton pair final states 
(e + e~, and r + r~) are very small. These contributions include the full 4-fermion 
final states, QCD-like final states and 77 scattering events. They are summed together 
in the table as a separate column. 

The SM describes the data well for each y/s and each step of the selection proce- 
dure, except for a single e/i candidate which was selected at y/s = 189 GeV. According 
to our Monte Carlo simulations only 0.019 events are expected from SM processes. 
This particular event is displayed in Figure |] for the r — <fi view and in Figure § for 
the r — 9 view. It has the following characteristics: 

• total visible energy E vis = 176 GeV (y/s = 189 GeV); 

• measured ECAL energy of the electron candidate: Eecal = (84.83 ±2.25) GeV; 

• measured momentum of the muon candidate: (84.5 ± 10) GeV; 
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Table 3: 95% confidence level upper limits on cr(e + e 
of a/s. 



e/x, er and /xr) as a function 



Channel 


e/i 


er 




Vs(GeV) 


<r[fb] 


<r[fb] 


<r[fb] 


189 


58 


95 


115 


192 < < 196 


62 


144 


116 


200 < v/s < 209 


22 


78 


64 



opening angle between the two tracks = 165°; 

one measured cluster in the forward detector with E = 7.2 ± 2.5 GeV. 



A kinematic fit with energy-momentum conservation is applied assuming that the 
forward detector cluster was due to a photon. The fit has a probability of 41%. 

Results of the present search are quantified in terms of 95% CL upper limits of the 
production cross section of LFV events, assuming uniform final state angular distri- 
butions. The upper limits are obtained with the method that takes into account the 
uncertainties on the signal efficiency and on the background expectation as explained 
in Reference [22|. These uncertainties include systematic contributions discussed in 
the next section. The upper limits are displayed in Table 3 for different intervals of 
a/s. Limits for different assumed angular distributions may be derived using the fact 
that the efficiency is uniform over the barrel region of the detector. 



5 Systematic Studies 

Five inputs are used to estimate the upper limits: the number of selected event 
candidates, the SM background contribution, the signal selection efficiency and the 
uncertainties on the background expectation and on the signal efficiency These un- 
certainties are estimated by repeating the analysis while varying the applied cuts, 
taking into account the following effects: 

• Efficiency: The two Monte Carlo methods (see Section 2) yield compatible 
results. The data-based method is applied only to the e/i channel, since at 
each y/s, the reconstructed r + r~ sample is too small to apply the method. A 
systematic error of 1% for e/i and of 2.1% for er and /ir is assigned based on 
the largest deviation between the methods. 

• Event selection: The following variations in the event selections are made, 
one at a time: the maximum number of allowed tracks per event is changed 
from 8 to 10; the minimum track momentum is increased from 250 to 500 MeV; 
the cone angle is opened from 15° to 30°; the cone minimum energy is increased 
from 10% to 20% of E beam ; the cut on the event thrust is set to 0.98 instead of 
0.95. 
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• Energy and momentum resolution: The electromagnetic calorimeter en- 
ergy resolution is about 3.5% for electrons with momenta above 70 GeV. The 
systematic error is estimated by scaling the measured electromagnetic energy 
of electron candidates by ±3.5%. The measured momenta of muon candidates 
are varied by ±10% to account for the momentum resolution. 

• Integrated luminosity: A 0.5% measurement error on the integrated luminos- 
ity is added to another 0.5% interpolation and averaging error between various 
grouped points(see Table 1), after which the corresponding error on the SM 
contribution is calculated. 

The different errors are added in quadrature. The final systematic uncertainty is 
about 3.5% on the efficiency and 5% on the background expectation. 



6 Conclusion 

We have no clear evidence for production of lepton flavour violating events such 
as e/i , er and /it production in e + e~ collisions between 189 GeV and 209 GeV. We 
observe one single e/i candidate, probably produced with initial state radiation, where, 
according to our Monte Carlo simulations, we expect 0.019 events to be produced from 
Standard Model processes. We have obtained the first upper limits for <r(e + e - — > 
e/i, er and /rr) as a function of y/s at LEP2 energies. The limits range from 22 fb to 
58 fb for the e/i channel, from 78 fb to 144 fb for the er channel and from 64 fb to 
166 fb for the e/i channel. 
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Figure 3: An r — view of the e/i candidate at 189 GeV. Two well measured tracks 
can be seen in the trackers (inner thick circle). The track in the right side deposited 
all its energy in the electromagnetic calorimeter (black trapezoid), and is the electron 
candidate. The track in the left side has a small electromagnetic energy deposit 
and has matching muon hits in the muon chambers (the arrow), and is the muon 
candidate. 
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Figure 4: An r — 9 view of the e/x candidate at 189 GeV. The beam direction is a 
horizontal line passing trough the intersection of the two tracks. The tracks are not 
back-to-back, and the missing momentum is compatible with the observed cluster in 
the forward detector which is the dark block close to the beam axis (right side). 



15 



